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ABSTRACT 

A  simplified  method  of  estimating  the  performance  of  >upersonic 
combustion  ramjet  engines  is  presented.  The  method  utilizes  stream 
thrust  concepts  and  enables  valid  periormance  estimates  to  be  made 
without  the  aid  of  a  computer  program;  only  a  few  simple  graphs  are 
required.  A  new  "reference''  stream  thrust  quantity  is  defined  and  showt. 
to  be  of  value  in  estimating  engine  flight  performance.  The  data  given 
in  this  report  enable  perfcnnance  estimates  .to  be  madr  for  hydrogen 
fueled  engines  operating  s to ichiometrica 11 y  for  speeds  in  excess  of 
about  Mach  8. 
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SYMBOLS 


Flow  Area 
Area  Ratio 
Reference  Area 
Speeo  of  Sound 
Inlet-Drag  Coefficient 

Characteristic  Velocity 

Specific  Heat  At  Constant  Pressure 

Internal  Thrust 

Universal  Gas  Constant  - 

Gravitiona!  Constant 

Static  Enthalpy 

See  Figure  19 

Total  Enthalpy 

Assigned  Enthalpy  of  Fuel 

Fuel  Specif’ c  Impulse 

Mechanical  Equivalent  of  Heat 

Ratio  of  Exit  Area  to  Inlet  Area  (A£/A0) 

Process  Efficiency 

Defined  in  Appendix  I 

Mass  Flow 

Fuel  Mass  Flow 

Air  Mass  Flo?' 
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SYMBOLS  (CONTD) 
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Subscripts 

Planes  0,  2,  {*).  3,  E,  R  are  shown  In  either  Figure  1  or  3. 
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SECTION  I 
INTRODUCTION 

This  report  Is  concerned  with  a  simple  approximate  analysis  of  the 
performance  of  scramjet  engines  operating  at  high  flight  speeds.  The 
method  of  analysis  was  conceived  In  1961  (Reference  1)  when  suitable 
computer  p!ograms  for  calculating  performance  were  not  available,  but  it 
was  not  published^  at  this  time  and  it  was  overtaken  by  the  development 
of  the  program  described  in  Reference  2.  Subsequently,  the  simplified 
method  became  of  interest  when  it  was  desired  to  study  the. effect  of 
perturbations  of  certain  component  efficiencies  on  engine  performance 
and  yet  avoid  a  major  production  run  on  the  computer.  This  method  of 
estimating  performance  Is  based  on  the  study  of  the  changes  in  stream 
thrust  of  the  engine  gas  flow  as  It  passes  through  the  engine. 
Consequently,  current  methods  of  defining  component  performance  are 
easily  introduced  into  the  analysis.  The  heart  of  the  method  lies  in 
the  fact  that  at  very  high  speeds  the  area  ratio  of  the  exhaust  nozzle 
is  large  and  the  state  condition  of  the  exit  flow  Is  normally  located 
in  the  non-dissociation  region  of  the  Molller  diagram.  This  fact  enables 
a  relatively  simple  approach  to  the  estimation  of  the  stream  thrust  of 
the  engine  gas  flow  at  the  nozzle  exit  to  be  made.  Using  this  simple 
approach,  engine  performance  at  speeds  of  about  Mach  15  have  been 
estimated  to  an  accuracy  of  one  to  two  percent  in  comparison  with  more 
sophisticated  computerized  calculations. 

This  report  is  submitted  in  the  hope  that  it  will  stimulate  the 
application  of  stream  thrust  concepts  to  the  analysis  of  engine  and 
component  performance. 
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This  report  is  a  more  general  version  of  Reference 
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EXPLANATION  OF  METHOD 


The  object  of  this  method  is  to  evaluate  the  specific  thrust  and 
fuel  specific  impulse  of  a  supersonic  combustion  ramjet  engine  operating 
at  a  given  flight  condition.  In  evaluating  these  performance  parameters, 
extensive  use  is  made  of  the  concept  of  air  specific  stream  thrust. 


Referring  to  Figure  1  the  internal  thrust  of  a  ramjet  engine, 
expressed  In  terms  of  the  stream  thrusts  of  the  Inlet  and  exit  gas 
streams,  is  given  by  the  conventional  expression: 


F  s  {stream -Thrust  j  -  {stream  Thrust  J  -  P0  (AE  -  A0) 


{■f-+PA 


L  -{-f*  paL  -  po  (ae  -  *0* 


The  specific  thrust  of  the  engine  may  be  written  as 


$a£  -  Sa0 


where  generally 


So 


J_  (M 
m0  l  9 


■f  PA 


} 


The  foregoing  expression  for  the  specific  thrust  of  an  engine  may  also 
be  written  as 


F 


(1) 


Ag 

If  attention  is  now  restricted  to  engines  in  which  the  ratio  j-  is  held 
constant  then  of  course  one  has  a 
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The  corresponding  fuel  specific  impulse  of  the  engine  is  defined  by 


On  inspection  of  Equations  fl)  and  (2),  it  is  apparent  that  for  a  given 
engine,  operating  at  a  specified  flight  condition  and  fuel/air  ratio,  the 
or: ly  unknown  quantity  is  Sag.  Once  this  quantity  is  determined,  the  engine 
performance  can  be  easily  evaluated.  The  unique  feature  of  this  simplified 
method  of  estimating  performance  is  that  it  permits  a  relatively  accurate 
determination  of  Sa^  without  resorting  to  complex  and  lengthy  calculations. 

Consider  an  engine  operating  at  a>given  condition  and  assume  for  the 
moment  that  the  flow  through  the  exit  nozzle  is  an  equilibrium  i?e.ntro"'tr 
process.  The  actual  expansion  through  the  nozzle  can  be  considered  as  a 
segment  of  a  hypothetical  isentropic  expansion  which  starts  from  the 
stagnation  condition  following  combustion,  and  proceeds  until  the  gas 
temperature  approaches  absolute  zero.  At  any  intermediate  state  condition, 
the  corresponding  air  specific  stream-thrust  can  be  calculated  from 

So  s  { I  4-  q } 

wnere  q,  T,  and  ft  are  known  «mi  V  can  be  obtained  from  the  Equation 

V2  =  2gJ  -  h) 

For  any  giver,  expansion  process  it  is  thus  possible  to  plot  the  vai iat^on 
in  Sa  as  a  function  of  a  suitable  reference  area  ratio  A/Ar  or  a 
reference  pressure  ratio  PyPg.  In  gas  dynamic  work,  it  is,  of  course, 
conventional  to  use  the  sonic  condition  as  a  suitable  reference  state. 
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Two  reference  values  of  specific  thrust  have  oftsn  been  used  In  ramjet 


performance  analysis,  namely  i.t*«rt*spf>nd i nn  t &  Uu*  sonic  condition  in 
the  expansion  and  Sa,^  which  is  thu  m  mum  spec; fig  stream  thrust 
obtained  by  expansion  to  an  atkt>‘>i!(le  **?> .?e»  <* ! »*•*  •«*'  -?»r5>.  rhese  y-iice1. 
are  indicated  on  figure  2  which  illustrates  a  tppi'gil  variation  of 
Sa/Sa^  with  A/A*.  In  the  performance  analysis  of  subsonic  combustion 
ramjet  engines,  this  type  of  chart  was  of  the  ut^m  value  This  was 
because  the  reference  area  A*  corresponded ‘is  the  geometrical  area  of 
the  choked  throat  of  the  exhaust  nozzle.  Thus  for  o  $fven  divergent  area 
ratio  Ap/A*..  it  was  possible  to  read  off  directly  the  sfc*utr$f  San -from 

the  chart.  However,  in  the  case  of  the  supersonic  comoustiofs  engine, 

;  *✓* 

this  tvne  of  chart  cannot  be  directly  utilized  since  the  sonic  reference 
area  Joes  not  correspond  to  any  specific  geometrical  area,  except  in  ths 
limiting  case  where  sufficient  heat  is  added  to  reduce  tne  exit  velocity, 


1 


*n 


4 


-* 

A 

I 


'S 

'fl 

i 


say  from  a  constant  area  combustor,  to  the. sonic  value. .  U  “-s,  frevertbs- 
less  informative  to  consider  the  application  of  this  chart -to  tbp-r 
scramjet  engine.  Referring  to'  Figure  3,  in  which  «v  acramjet  with,  a.-; 
constant  area  combustor  is  considered,  it.  is  apparent  that  the  specific 
stream  thrust  at  plane  3  is  determined  by  the  flow  conditions  at  thsr 
diffuser  discharge,  the  stream  thrust  contributed  by  the  fuel  injection 
system,  ano  the  frictional  losses  in  the  combustor.  If,  for  the  time 
being,  the  Utter  two  effects  are  neglected,  we  have 


So 


3 


Sa2 


iyv 

~vt 


and  thus  for  this  particular  example  the  s?ecific  stream  thrust  after  -3 

t 

j 

i, 
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Figure  2.  Variation  sf  Impulse  Ratio  with  Area  Ratio 
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Figure  3.  Hypothetical  Expansion  Process 
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combustion  is  uniquely  determined  once  the  operating  condition  of  the 

intake  is  known.  Th^s  for  a  specific  flight  condition  and  adiabatic 

inlet  flow,  Sa-j  is  dependent  only  on  the  amount  of  diffusion, 

performed  by  the  intake  and  Is  independent  of  the  efficiency  of  the 

intake  (real  gas  and  viscous  effects  being  neglected).  It  must, 

however,  be  noted  that,  for  a  fixed  value  of  Ar/A  .  the  nozzle  area 

An  Ac  ^  ^ 

ratio  Is  given  by  *  r-  .  Tp,  thus  for  a  given  amount  cf  diffusion 

4  *2  o 

the  effect  of  reducing  inlet  efficiency  is  to  decrease  ine  available 
area  rctfo  of  the  nozzle. 


At  flight  speeds  below  about  Mach  7.5,  it  is  i^ssible  to  utilize 
charts  of  the  type  i  i lustrated ,i.i  Figure  2  to  determine  the  performance 
cf  scramjet  engines.  In  this  case  the  chart  is  entered  at  the  value  of 
Sa3/Safc,  where  Sa3  corresponds  to  the  specified  inlet  condition,  and  a 
value  of  Aj/A^  is  obtained.  The  chart  is  then  re-entered  at  an  area  ratio 

AE  9  As 

A«  As  A» 

and  the  corresponding  value  of  Sa^  obtained.  By  ris  procedure  it  is 
possifc-le  to  utilize  cherts  wh»ch  are  already  available  for  subsonic  com¬ 
bustion  engines  for  calculating  the  performance  of  scramjets.  However, 
such  charts  ars  nGt  readily  available  for  speeds  in  excess  of  about 
Mach  *\fi  and  considerable  labor  is  required  to  produce  the  necessary 
charts  to  cover  all  probable  conditions  of  engine  operation.  Further¬ 
more,  as  higher  flight  soeeds  are  considered,  it  becomes  increasingly 
difficult  to  calculate  A#  and  P*  since  the  sonic  state  points  are  located 
in  the  high-enthalpy  regions  closely  corresponding  to  the  stagnation 
condition  cf  the  flow.  Indeed  a  cor <ii tier,  is  soon  reached  when  the 
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sonic  state  poiits  correspond  to  enthalpies  which  exceed  the  range  of 
available  Mol  her  diagrams.  Additionally  this  procedure  for  determining 
performance  is  physically  unattractive  since  the  flow  within  the  engine 
never  approach'  s  the  sonic  reference  condition.  In  fact,  the  state 
points  corresponding  to  the  actual  flow  through  the  engine  are  usually 
located  in  the  lower  regions  of  the  Moilier  diagram.  It  is  also  relevant 
to  note  that  at  very  high  speeds  the  area  ratio  of  the  exhaust  nozzle 
becomes  very  urge  and  the  state  condition  of  an  exit  flow  which  results 
from  an  equii  brium  expansion  is  normally  located  in  the  non-dissociation 
region  of  the  Moilier  diagram.  This  latter  consideration  permits  £he 
use  of  a  simple  approximate  approach  to  the  estimation  of  Sa^* 

The  basis  of  this  new  method  is  to  use  a  reference  specific  stream 
thrust  SaR  which  in  concept  corresponds  more  closely  to  Sa,,^  rather  than 
Sa#;  one  then  works  back  upstream  to  determine  Sag  rather  than  downstream 
from  the  value  of  Sa*  as  previously  discussed.  Of  course,  Sa^*  itself 
cannot  be  usea  as  a  reference  condition  since  it  corresponds  to  an 
infinite  expansion  area  rotio. 

In  this  report  the  reference  specific  stream  thrust  arbitrarily 
defined  as  the  value  corresponding  to  an  expansion  to  a  static  tempera- 
cure  of  39Q°R  (216.5eK);  this  te.c<perature  corresponds,  of  course,  to  the 
standard  tropooause  value  and  is  a  rational  ’sink"  temperature.  The 
reference  specific  stream  thrust  is  easily  calculated;  thus  for  an 
adiabatic  equilibrium  expansion  from  post-combustun  conditions 

2  i.  2 

V*  Vd 

M.  =  h,  4  c  h-  4  “ — 

♦s  3  StgJ  «  2gu 

so  tnat 
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Where 


and  thus 


V  =  2*J  ■  (Htj  -  ».„) 


1  \  frh} H?  +  dr  K + 

s.  ,  (l+„  {-f  ♦ 


Tie  value  of  Samax  may  be  calculated  in  a  similar  fashion.  Computed 
values  of  Sa^,  expressed  for  accuracy  as  -  V0/g,  for* stoichiometric 

hydrogen/air  combustion,  are  given  in  Figure  4  for  various  combustion 
efficiencies.  Additionally,  values  of  Samax  and  Sa#  are  compared  in 
Figure  0  for  stoichiometric  products  and  r,c  *  1.00. 

For  any  given  stagnation  condition  it  is  a  simple  matter  of  calcula¬ 
tion  to  obtain  the  variation  of  Sa  with  area  for  a  simple  adiabatic 
equilibrium  expansion.  Thus  at  any  state  point  in  the  expansion,  P,  T, 

r-  1 

and  R  =  ^  are  known,  and  the  corresponding  flow  velocity  can  be  calcu¬ 
lated  from  the  enthalpy  change.  Thus 

So  =  (I  +  0)  (t  +  JH 

and  the  flow  area  per  unit  flow  can  be  calculated  from 

A  .  RT 
(I  ♦  q )  m0 

It  is  thus  possible  to  produce  curves  of  Sa/Sap  versus  A/AR  for  any 
initial  stagnation  condition  and  these  curves  are  simi  to  those  of 
Si/Sa*  versus  A/Ar  illustrated  in  Figure  2.  For  *  given  ^uivaience 
ratio  the  Sa/Sa^  '  A/Ag  curves  are  only  dependent  t  *  i’e  total  enthalpy 
of  the  inlet  stream  and  the  initial  pressure;  however,  as  the  expanjicn 
proceeds  away  from  stagnation  conditions  the  dependence  of  these  curves 
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on  the  absolute  pressure  level  disappears.  Graphs  showing  the  actual 
variation  of  Sa/Sap  with  A/Ar  are  shown  in  Figures  i  to  18;  these  graphs 
are  for  the  stoichiometric  combustion  of  hydrogen  air  mixtures  at  initial 
conditions  of  1  atm.  pressure.  Superimposed  on  these  curves  are  values 
of  the  corresponding  expansion  pressure  ratios. 

The  above  curves  can  be  utilized  to  calculate  the  specific  stream 
thrust  at  tne  exit  of  a  scramjet  engine  nozzle.  For  example,  consider  an 
engine  with  i  simple  constant  area  combustor  operating  at  an  equivalence 
ratio  of  l .0. 

¥ 

From  Figure  4,  Sa^  miy  be  obtained  from  the  value  of  Sa^  -  ~  corre¬ 
sponding  to  the  flight  condition.  The  quantity  $a3  will  be  known  for  a 
given  intake  operating  point  (i.e.,  Vj/Vq  giver;  and  correspondingly  Ag/A 
will  be  r-  wn  for  a  given  engine  geometry,  thus  A£/A3  =  Ac/Aq  .  Aq/Aj. 

The  c  .it  stream  thrust  may  now  be  found;  thus  entering  the  charts  in 
Figures  £  >.tc. ,  at  a  value  of  Sa/Sap  *  Saj/Sap  corresponding  to  the  flight 
condi tioc  one  cbtains .A3/A3;  tnis  latter  value  may  be  used  to  calculate 
Ar/AR  *-  : -/A3  .  A3/Ar.  The  chart  is  now  re-entered  at  A/Ap  =  Ajr/Ap  and 
Sop/Sag  wn  be  read  off,  and  Sac  found. 

The' above  process  is  extremely  simple  when  constant  area  combustors  are 

consider  d  and  it  is  relatively  easy  to  study  the  effects  of  varying  sucn 

parameters  as  inlet  efficiency,  fuel  injection  conditions,  and  combustor 

friction  losses  on  engine  performance.  Non-equilibrium  and  other  nozzle 

Ics'is  can  also  be  introduced  into  the  analysis  provided  such  losses  are 

expressed  in  terms  of  stream  thrust  changes.  The  method  is  very  accurate 

for  engines  operating  at  high  flight  speeds  {M  >  12}  with  correspondingly 

high  expansion  ~>tios  in  the  nozzle.  However,  before  discussing  the  validity 
and  »i  .  tat  ions  this  method  of  approach,  it  is  considered  desirable  to 

i’lusfcetp  the  nethod  further  by  discussing  sample  performance  calculations , 
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SECTION  III 

PERFORMANCE  CALCULATIONS 

The  ~ey  item  to  be  calculated  In  assessing  performance  Is  the  value 
of  Sa^  which  in  turn  Is  dependent  on  Sa^  and  the  available  expansion 
area  ratio  A^/Aj.  To  determine  these  latter  two  quantities  it  is 
necessary  to  calculate  tne  Intake  area  ratio  and  the  variation  of  Sa 
through  the  intake  and  combustor. 

1.  Calculation  of  Intake  Area  Ratio  and  Exit  Specific  Stream  Thrust 

As  noted  in  the  Introduction,  for  a  fixed  value  of  Ag/AQ  the 
area  ratio  of  the  nozzle  is  dependent  on  the  area  ratio  of  the  intake. 
The  area  ratio  of  the  intake  is  easily  calculated  once  the  degree  of 
diffusion  and  the  associated  efficiency  of  diffusion  are  known.  For  a 
simple  adibatic  compression  process  from  velocity  VQ  to  velocity  V2 

"o  +  T&r  *  **  ♦  “S3" 

So  that 


For  assumed  values  of  yr-,  h2  is  easily  calculated  for  given  flight 

o 

condition  using  the  tables  in  Reference  3.  The  Intake  efficiency  may 
be  defined  in  various  ways  (See  Reference  4)  but  if  the  kinetic  energy 
definition  is  chosen  then  from  Figure  19,  one  has 


tNote  that  conventionally  nKE  is  often  expressed  as 

\E  =  S  +  <<-  S' 
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Where 


•Ht  =  ho  +  -ip- 


So  that 


hO  s  “o  +  d 


KE' 


The  area  ratio  corresponding  to  the  compression  process  is 


*o 

A* 


where  *=•  is  the  only  unknown.  Along  the  isentrope  corresponding  to 

0  p2 
$2  (Figure  19)  the  ratio  of  the  pressures  -p= — ppy  can  be  found  from 

the  tables  of  Reference  3  using  the  relative  pressure  function.  These 
tables  will  normally  cover  all  practical  conditions  of  intake  operation. 


The  value  of  Sa2  can  be  calculated  from 


-  ■ 


2.'  Variation  of  Stream  Thrush  Through  Combustor 

For  the  present,  attention  is  restricted  to  the  constant  area 
combustor;  thus  any  change  in  the  stream  thrust  will  arise  due  to  fuel 
injection  or  frictional  effects.  Considering  first  fuel  injection, 
possible  geometries  are  sketched  in  Figure  20.  If  wall  slot  injectors 
are  used  then  the  area  ratio  available  for  expansion  in  the  nozzle  is 
reduced,  but  of  course  there  is  a  contribution  to  the  stream  thrust  at 
exit  from  the  combustor.  For  fully  axial  injection  of  fuel  at  velocity 
Vf  and  temperature  T^,  the  corresponding  fuel  stream  thrust  is 
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and  the  corresponding  air  specific  stream  thrust  of  the  injected  fuel  is 


to. 


rv«  ,p*An .  „  /v*  *  RfTf 

nT+*n ' ’vr*— 


In  Figure  21  the  variation  of  V f  and  Tf  for  various  fuel  total  temperatures 

r  ~i 

is  shown  and  in  Figure  22  the  variation  in  Saj^  with  Vf  is  given.  The 

*>  —  -  J 

area  required  for  fuel  injection  can  be  calculated  from  the  expression  •< 


R2  T,  Pf 


pf  =  p2 


1= , 


R,  T(  v2 


A2  R2  T2  Vf 

The  specific  stream  thrust  at  Plane  3  is  thus  given  by 


So,  -  So»  + 


Whs 


For  the  strut  injection  configuration  shown  in  Figure  20  there  is  no 
increase  in  the  combustor  area  from  Plane  2  to  Plane  3.  However,  it  is 
necessary  to  account  for  the  drag  force  on  the  injector  when  calculating 
the  exit  stream  thrust. 

Although  friction  losses  can  be  appreciable,  no  attempt  will  be  made 
here  to  give  a  generalized  method  of  calculating  such  losses  since  the 
computation  of  friction  losses  depends  significantly  on  the  configuration 
of  the  combustor- injector  combination  and  on  the  precise  nature  of  the 
boundary  layer  flow.  However,  friction  force  losses  can,  of  course,  be 
treated  in  a  parametric  study. 

The  extension  of  this  analysis  to  include  the  constant  pressure 
combustor  is  given  in  Appendix  I. 
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Vf  *  10" 5  ft/*« 


Figure  22.  Air  Specific  Impulse  for  Hydrogen  Fuel  vs  Fuel  Velocity 
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3.  Calculation  of  IM  Perfof'r^nce  of  a  Non-Ideal  Nozzle 
The  value  of  Sag  given  in  the* charts  in  figures  6  to  18 
corresponds  to  ao  iser.tropic  equilibrium  expansion  through  a  given  area 
ratio.  The  performance  of  a  real  nozzle  is  reduced  relative  to  the  ideal 
nozzle  because  of  non-equilibrium  effects,  friction  losses,  divergence 
losses,  and  other  sources  of  entropy  production.  Usually  such  losses  are 
accounted  for  by  coefficients  applied  to  the  ideal  results.  A  large  num¬ 
ber  of  coefficients  have,  been  defined  for  nc-zzle  analysis  but  for  the 
purposes  of  this  paper  only  these  based  on,  or  directly  related  to,  the 
thrust  produced  by  the  exit  gas  stream  are  considered.  These  coefficients 
are  very  easily  introduced  into  the  analysis.  For  example,  one  common 
coefficient  is  defined  as: 


=  The  ratio  of  the  actual  nozzle  vacuum  thrust  to  the 
ideal  vacuum  thrust  of  an  isentropic  nozzle  of  the 
same  area  ratio,  the  mass  flow  and  applied  pressure- 
ratio  being  the  same  in  both  cases. 

Thus 

/Pc)  rSftp) 

_  v_  cctuol  _  1  c,octuol 

*  -  IdVoi 

Another  coefficient  frequently  encountered  is  similar  to  ->s  but  it  is 
based  on  the  nozzle  gauge  thrust.  Thus: 


=  The  ratio  of  the  gauge  thrust  of  the  actual  nozzle  to 
the  gauge  thrust  of  an  ideal  nozzle  of  the  same  area 
ratio,  tne  macs  flow  and  applied  pressure  ratio  being 
the  same  in  both  cases. 
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Some  sample  calculations  were  made  to  compare  the  accuracy  of  this 
approximate  method  with  the  more  detailed  computer  program  calculations 
given  in  Reference  2,  which  were  basically  one-dimensional,  equilibrium  flow 
calculations. 

The  results  are  shown  in  Figure  23.  Four  sets  of  calculations 
were  made  for  speeds  of  12,000,  14,000,  16,000,  and  18,000  ft/sec  for  a 
stoichiometric,  hydrogen-fueled  engine.  An  intake  performance  of  Kq  *  0.95 
was  assumed  and  an  isentropic  nozzle  expansion  with  ns  =  1.0.  The  altitude 
corresponded  to  a  dynamic  pressure  of  500  Ib/sq  ft. 

The  agreement  between  the  approximate  and  the  inore  detailed 
calculations  is  seen  to  be  very  satisfactory  and  the  agreement  appears  to 
improve  with  increasing  flight  speed. 

A  lower  speed  check  case  was  also  calculated  for  a  Mach  9  engine 
with  =  0.975,  -:c  ~  0.95,  and  A^/Aq  a  1.5.  The  results  are  shown  in 
Figure  24  and  once  again  the  agreement  between  the  approximate  and  the 
detailed  method  is  good. 
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Velocity  Ratio  V2/V0 

Figure  24.  Comparision  of  Exact  and  Simplified  Calculations 
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A  specific  example  of  a  simple  slide  rule  calculation  which  is 
intended  to  illustrate  the  use  of  the  performance  curves  is  given  in 
Appendix  II. 
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SECTION  IV 

INCORPORATION  AND  ANALYSIS  OF  EXPERIMENTAL  RESULTS 

Another  use  of  the  analysis  presented  earlier  Is  the  interpretation 
of  the  results  of  component  tests  in  terms  of  overall  engine  performance, 
for  example,  if  an  experimental  combustor  is  considered,  then  it  may  be 
desirable  to  "fit"  a  hypothetical  intake  and  nozzle  to  the  combustor  to 
determine  the  effect  of  combustor  performance  variation  on  engine  perform¬ 
ance  at  a  given  flight  condition. 

In  assessing  the  performance  of  components  it  is  evident  that  the 
important  qualities  to  establish  are  the  actual  thrust/drag  forces  on  the 
components.  Although  the  basic  thermodynamic  efficiencies  of  components 
were  initially  of  interest  in  development  and  in  cycle  analysis,  the  trend 
of  assessment  has  been  toward  direct  measurement  of  the  forces  acting  on 
components.  Thus  the  experience  gained  in  previous  years,  both  in  the 
testing  of  ramjet  combustors  and  of  freejet  rigs,  has  been  applied 
to  the  measurement  of  scramjet  inlet  performance  (see  References  5  at' 

6).  Similarly,  one  may  expect  such  techniques  to  be  applied  to  both  com¬ 
bustor  and  combustor- nozzle  development.  Apart  from  the  desirability  of 
direct  force  measurements,  this  trend  has  been  forced  by  the  inadequacy 
of,  and  anomalies  associated  with,  J  -  ;se  of  point  located  instruments 
to  assess  average  flow  properties. 

It  will  th-  j  be  apparent  that  the  determination  of  forces  acting  on 
components,  by  means  of  force-balance  test -rigs  alliec  with  mass  flow 
measurements,  permits  direct  calculation  of  me  changes  in  air  specific 
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stream  thrust  associated  with  a  component.  The  main  point  to  be  made 
here  is  that  whether  the  stream  thrust  change  is  estimated  theoretically 
or  measured  experimentally ,  it  is  a  relatively  easy  matter  to  introduce 
specific  component  performance  Into  this  analysis  or  alternatively  to 
use  the  analysis  as  a  means  to  express  the  impact  of  component 
performance  on  flight  performance. 
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SECTION  V 
CONCLUSION 

A  simple  method  of  calculating  the  performance  of  supersonic 
combustion  engines  operating  at  high  flight  speeds  has  been  devised. 

The  method  has  been  shown  to  yield  acceptable  results  for  speeds  in 
excess  of  9,000  fps  but  no  attempt  has  been  made  to  assess  the  validity 
at  lower  flight  speeds.  Charts  are  only  given  for  stoichiometric 
hydrogen/air  combustion  since  hydrogen  fuel  is  the  prime  candidate  for 
high  speed  engines  and  operation  at  lower  equivalence  ratios  is  not 
desirable  at  high  flight  speeds. 

The  concept  of  a  "reference"  station  corresponding  to  adiabatic 
expansion  from  post  combustion  stagnation  conditions  to  a  reference 
temperature  condition  may  find  further  application  in  other  engine  cycle 
studies. 
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APPENDIX  I 

STREAM  THRUST  ANALYSIS  OF  CONSTANT  PRESSURE  COMBUSTOR 

The  calculation  of  the  stream  thrust  at  exit  from  a  constant  pressure 
combustor  is  relatively  straightforward  if  fuel  ‘injection  momentum  and 
skin  friction  effects  are  neglected.  Using  these  assumptions,  the  method 
of  approach  is  as  follows.  One  can  write  for  the  constant  pressure 
combustor 


For  given  inlet  conditions,  the  only  unknown  in  the  above  equation  is  the 
combustor  area  ratio  A3/A2.  This  ratio  can  be  expressed  from  the  con¬ 
tinuity,  momentum,  and  energy  equations.  Continuity  yields 

(I  q)p2A2V2  =  />3 A3V3 


3  -  1 

s - \ 


The  momentum  aquation  yields 


0  ♦  q)  V5  *  V2 


=  <.♦ 


From  the  energy  equation  one  has 


h3  =  Hf_  - 


3  2gJ 
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where 


s  ■  {tttH  +  i^ 

Thus  for  given  assumptions,  is  easily  computed  as  a  function  of  fuel 
conditions  and  flight  speed  and  h3  can  be  computed  for  given  combustor 
entry  conditions  by  using  Equation  1-2.  A  relationship  between  hq  and 
Ty/^  is  shown  in  Figure  25.  (This  relationship  is  exact  for  a  pressure 
of  1  atmosphere;  minor  error  is  involved  for  use  at  other  pressures.) 

This  graph  can  now  be  used  to  calculate  A3/A2.  Finally,  Sa3  can  be 
obtained  from  Equation  1-1. 

The  above  approach  can  be  extended  to  cases  where  the  axial  fuel 
momentum  term  is  significant.  To  speed  up  such  calculations,  the  follow¬ 
ing  analysis  is  used  to  derive  some  useful  charts.  Again  let  us  utilize 


the  energy  equation 


V* 

h3  s  H»,  “W 


whence 


TrqH?  +  h0  ♦  V02 /2«J 


+  h° 

H»  =  — - + - s  - 

*3  I  +  q  »  +  q  I  ♦ 


For  an  assumed  fuel  temperature  of  20.39°K  and  stoichiometric  combustion, 
one  has,  where  all  enthalpies  are  now  based  on  Reference  7, 

Hj  =  60,309  BTU/lb. 
q  *.  0.02928  lb-fuel/ib-air 


T0  =  390 °R  ,  h0  =  203.6 


then,  assuming  nc  =1.00 


1913.4  ♦ 


2gJ(  I  +  <t ) 
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V0 _ _  VS 


From  the  momentum  equation 

V2  +  q^V2  =  (l+q)V, 


so  that 


v3  ‘ 


one  then  obtains 


V, 

V?  ♦  <1  ^  v2 


*  ,9*3  4  +  2flJ  n  ♦  q)  “  I  +  q 

ZqJ 


U  *  >913.4  +  ■■"  -r  I 

3  2gJ  (I  +  q) 


Writing 


one  has 


(•••?} 


i  +  q 
I  +  q 


Q"^r) 


J  ♦  q 


h3  *  ,913.4  +  K^p 

Figure  26  shows  the  variation  of  hj  with  V2/V0  for  various  free  stream 
velocities  and  constant  pressure  combustion  with  no  fuel  injection.  Figure 
27  is  a  plot  of  K  versus  V^/Vg  for  various  values  of  fuel  injection  velocity 
ratio.  Figure  25  has  already  shown  the  variation  of  versus  h3- 
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Figure  26.  Combustor  Enthalpy  -  Inlet  Velocity  Ratio 


47 


■  —  -  ->*•  t*  - j. Ate*.  9&**^^*?***  1  ^rr>fifi4rri  !*»■■  “ 


Figure  27.  factor  for  Computing  Combustor  Enthalpy  vs 
Inlet  Velocity  Ratio 
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The  next  quantity  required  is  the  area  ratin'  of  the  combustor.  Using 
the  continuity  equation  we  have 


A-  fgt 

-t*  *  ( I  ♦  q) - 

*2 


V* 


and  making  use  of  the  equation  of  state 


«  (I  ♦«)" 


f  h  \ 

T*  i 

VP.^sV 

,  one  has 

h 

*  « 

\ 

% 

vf 

1 

V*  . 

which  is  the  required  expression. 

One  must  now  determine  the  specific  stream  thrust  at  station  3. 
This  can  be  accomplished  by  once  again  making  use  of  the  momentum 
equation  and  the  definition  of  air  specific  stream  thrust. 

The  momentum  equation  applied  to  the  combustor  yields 


<*,»,  +  P,A,-  Vs  *  + 


Now  from  the  definition  of  air  specific  stream  thrust,  and  the  fact  that 
pi  sps!pf 

P, 

+ 

’2  mO 


So,  s  So2  +  [So]H  +  -^  {a,-  (A2  +  Af)} 

3  * 

and  utilizing  the  continuity  equation 


s«s  •  So,  >  (:  +  £)} 
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finally 


Thus  with  a  knowledge  of  AyAg,  it  is  possible  to  calculate  Sa^  and 
then  to  proceed  to  the  calculation  of  iaE.  Finally,  the  overall  engine 
performance  can  be  determined. 
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APPENDIX  II 

EXAMPLE  OF  THE  USE  OF  PERFORMANCE  CURVES 


Consider  a  hydrogen  fueled  engine  operating  with  stoichiometric 
fuel/air  ratio,  n.  *  1.0,  at  a  flight  speed  of  16,000  fps  and  ambient 
temperature  500°R.  Assume  nKE  *  0.98,  r,N  *  1.00,  and  fuel  Injection 
normal  to  main  airstream.  A  constant  area  combustor  is  assumed  (A2  *  A3) 
and  also  we  assume  AQ  *  A^.  We  shall  calculate  the  engine  performance 
fcr  an  intake  velocity  ratio  V.>/V0  c  0.95. 

Initial  Conditions 


Given 


V0  16,000  fpt  Tq  *  500*R 


=  497.3 


=  5112.5 


V0  RTb 

So-  s  -  + - -  r  499 

0  g  V- 


From  Figure  5, 


Vo 

”  *5“  *  76.6  , 


Sag  «  573.9 


From  Figure  6, 


S°ma*  o’  *  ,  $°?r.aK  *  577,5 


Intake  Analysis 


ho  *0  *  (*  ‘  VKt)  zgj 


From  Reference  3 


x'o  *  119.5  +  002  -^3- 


*  221.75 
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From  Figure  14, 


Assuming 


aj/aR  9  7.0  x  10“5 


ae  a0 

at  x -at  8  7  96 


From  Figure  13, 


whence 


a£  _  *o  As  .4 

*»  ■  t;  Ti  * 5  57  - 10 


SOe 

fsj  *  a»* 


SOe  «  523.4 


Performance 

For  the  engine  geometry  assumed,  where  A-  »  Art,  we  have  from 
Equation  (2)  *  0 

X*P  *  "q*(SoE  ~  $*0) 

«  2A.4 

6.02928 

s  834.5  tec». 
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